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SECTION  I 


INTRODUCTION 

This  report  uses  the  modified  air  and  water  filtration  theory 
from  CEEDA-TR-77-1  to  develop  a common-based  rational  design 
sequence  for  graded  media  water  filters.  Characterization  of  these 
systems  must  include  terms  for  all  mechanisms  having  significant 
influence  upon  the  relationships  between  flow,  time,  pressure,  and 
efficiency. 

Since  the  model  developed  is  new,  its  application  to  real 
filtration  systems  is  unknown.  Furthermore,  the  sensitivity  of 
the  design  variables  that  characterize  the  suspended  particles  in 
the  water  is  unknown  because  these  data  are  not  used  with  current 
water  filtration  theories.  For  these  reasons  a subsequent  report 
will  evaluate  the  model  for  data  sensitivity  and  accuracy  from 
field  data. 


SECTION  II 


DISCUSSION 

The  model  requires  information  concerning  the  sand  (or  other 
graded  media)  characteristics,  the  fluid/particle  suspension,  and 
the  filter  operating  conditions.  Using  these  data,  the  model 
determines  the  statistical  characteristics  of  the  sand  bed  and  the 
suspended  particles,  the  effective  particle  size  distribution  through 
the  sand  bed,  the  single  sand  grain  collection  efficiency,  the  single 
suspended  particle  collection  efficiency,  and  the  sand  bed  solidarity 
factor.  At  this  point  the  model  integrates  numerically  to  determine 
the  pressure  drop  caused  by  the  sand,  the  suspended  particles  col- 
lected in  the  sand,  and  the  filter  efficiency  for  both  the  sand  and 
the  suspended  particles  (the  cake)  at  given  time  increments.  In 
addition,  it  computes  the  average  efficiency  from  the  start  to  any 
point  in  the  filter  cycle.  The  characteristics  of  the  sand  bed  and 
the  filter  cake  are  then  modified  for  each  increment  to  reflect  the 
net  change  in  their  characteristics  caused  by  particles  collected 
during  the  previous  integration  increment.  Figure  1 is  a logic  flow 
diagram  of  the  compute'  program  used  to  accomplish  the  above  described 
mathematical  sequer  he  logic  and  mathematical  expression  contained 
in  Figure  1 are  di  herein. 

To  facilitat  ter  understanding  of  the  model,  input  informa- 

tion concerning  , icer  operation  includes  flow  rate  per  unit  area  of 
filter,  total  filter  cycle  time,  integration  increment  cycle  time 
(accuracy  of  the  numerical  integration  technique  is  inversely  propor- 
tional to  the  integration  increment),  fluid  temperature,  fluid  den- 
sity, and  fluid  dynamic  viscosity.  The  design  variables  requiring 
definition  for  the  sand  bed  include  mass-mean  sand  size,  log-normal 
standard  deviation  of  sand  size,  filter  sand  bed  depth,  density  of 
the  discrete  sand  grains,  average  bulk  density  of  the  sand  bed,  and 
maximum  bulk  density  of  the  sand  bed. 

For  a specific  sand  the  mass-mean  size,  bulk  density,  maximum 
bulk  density,  and  discrete  sand  grain  density  must  be  evaluated  in 
the  laboratory.  These  measurements  are  performed  using  sieve  analysis 
techniques.  The  sand  size  is  plotted  in  percent  by  mass  passing  the 
sieve  as  a function  of  sieve  size.  From  these  data  the  remaining 
particle  size  distribution  variables  may  be  calculated  from  the  math- 
ematical characteristics  of  a log-normal  distribution.  A more  complete 
discussion  of  sieve  analysis  for  filter  sands  and  determination  of  sand 
bed  characteristics,  such  as  bulk  density  and  discrete  sand  grain 
density,  has  been  written  by  several  investigators  (References  1-7). 

Determination  of  the  characteristics  of  the  suspended  particle 
population  is  also  necessary  for  the  filtration  model.  The  method 
used  to  accomplish  this  may  vary  according  to  the  particle  size  range 
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Figure  1.  Logic  Flow  Diagram  for  the  Proposed 
Sand  Filtration  Model 
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INPUT  SAND  MEDIA  CHARACTERISTICS: 

Mass  mean  sand  size  (cm) 

Log-normal  standard  deviation 
Filter  depth  (cm) 

Discrete  sand  grain  density  (g/cm3) 
Average  sand  bed  bulk  density  (g/cm3) 
Maximum  sand  bed  bulk  density  (g/cm3) 


CALCULATE  SAND  BED  GROSS  CHARACTERISTICS: 


1 . Bed  depth  increment  at  10%  total  depth  (cm) 

2.  Incremental  bulk  density  change  per  layer 

3.  Top  layer  bulk  density  - maximum  (g/cm3) 


IFB  > 10 


DO  LOOP  15 
IFB  = 1 to  10 


IFB  < 10 


CALCULATE  STATISTICAL  CHARACTERISTICS  FOR  SAND  BED: 

1 . Geometric  mean  (cm) 

2.  Average  mean  (cm) 

3.  Surface  mean  (cm) 

4.  Effective  size  as  a collector  (cm) 


Calculate  suspended  particle  diffusion  coefficient,  D. 


Figure  1.  (Continued) 
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DO  LOOP  16: 
IFB  = 1 to  10 


IFB  > 10 


CALCULATE  CAKE  CHARACTERISTICS: 


. Friedlander  efficiency 
. Solidarity  factor 


CALCULATE  CONSTANTS  FOR  SAND  LAYER  IFB: 


1.  Sand  grain  effective  size  (cm 

2.  Sand  grain  drag  coefficient,  K2S 

3.  Single  grain  collection  efficiency 

4.  Sand  bed  solidarity  factor 


OUTPUT  BASIC  FILTER  DATA 


AND  INPUT  DATA 


DO  LOOP  10: 

ITIM  = ITIM2,  ITIM1 , ITIM2 

i .e. , start  0 ITIM2 

increment  0 ITIM2 
until  0 ITIM1 


ITIM  < ITIM1 


INCREMENTAL  INTEGRATION  during  filter  run 
for  pressure,  flow,  time,  and  efficiency 


CALCULATE  NEW  VALUES  FOR  PARTICLE  CAKE: 
TT  Calculate  new  cake  thickness 

2.  Calculate  new  cake  solidarity 

3.  Pressure  drop 

4.  Efficiency 

5.  Reduction  in  suspended  solids  load 


Figure  1 (pontinued) 


2.  Pressure  drop  across  suspended  solids  cake 

3.  Pressure  drop  across  sand  bed 

4.  Total  pressure  drop 

5.  Efficiency  across  suspended  solids  cake 

6.  Efficiency  across  sand  bed 

7.  Total  efficiency  at  time  t 

8.  Average  efficiency  at  time  t into  filtration  run. 


Figure  1 (Concluded) 
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contained  in  the  population.  A primary  method  for  particle  size  deter- 
mination could  be  hydrometric  analysis.  This  method  works  well  for 
particles  as  small  as  1.0  urn  in  diameter  (Reference  8).  Secondary 
methods  such  as  zonal  centrifugation  and  nuclei  counting  techniques 
(Reference  8)  may  be  employed  if  a substantial  mass  of  the  particle 
population  is  below  1.0  ym  in  diameter.  Suspended  particle  bed  bulk 
density  and  discrete  particle  density  can  be  determined  from  gravi- 
metric analysis  (Reference  9). 

With  this  information  the  sand  grain  or  the  suspended  particle 
size  distribution  can  be  classified.  These  calculations  include 
determination  of  geometric  mean  particle  size: 

dg  = exp[ln  dj^  - 6.908  In2  (og)]  (1) 

where  dq  is  the  geometric  mean  particle  size,  d^  is  the  mass  particle 
size,  and  og  is  the  log-normal  standard  deviation.  Next,  the  surface 
mean  particle  size  is  determined: 

dsu  = exp[ln  d^  - 1.151  ln2(og)]  (2) 

where  dsu  is  the  surface  mean  particle  size.  The  arithmetic  mean 
particle  size  is  then  calculated: 

dA  = exp[ln  d^  - 1.263  In  (og) ] (3) 

where  dA  is  the  arithmetic  mean  particle  size.  When  this  is  accom- 
plished, the  effective  particle  size  as  a collector  is  calculated  from 
Chen's  (Reference  10)  equation: 


where  de  is  the  effective  particle  size  for  the  collector. 

At  this  point,  all  of  the  basic  input  data  have  been  acquired 
and  the  basic  characteristics  of  both  suspended  particles  and  the 
sand  bed  have  been  calculated  (see  Figure  1).  The  first  filtration 
design  variable  to  be  calculated  is  the  particle  drag  coefficient  as 
follows: 


where  K2p  is  the  particle  drag  coefficient  (sec/cm),  pBp  is  suspended 
particleKbulk  density  (g/cm3),  d is  the  effective  particle  size 
(cm),  and  Pp  is  the  discrete  particle  density  (g/cm3).  The  next  cal- 
culation determines  the  suspended  particle  diffusion  coefficient.  The 
geometric  mean  suspended  particle  size  is  used  for  this  calculation  as 
suggested  by  Chen  (Reference  10)  and  Fuchs  (Reference  11): 


n - 1.2096  x 10‘17  T 

' — Tt <6 

g gp 

where  D is  the  diffusion  coefficient  (cm2/sec),  T is  temperature 
(degree  Kelvin),  yg  is  dynamic  viscosity  (poise  = g/sec-cm),  and  d 
is  the  geometric  mean  suspended  particle  size  (cm).  9P 

Subsequent  to  these  calculations,  the  sand  bed  is  divided  into 
ten  equal  depth  increments.  The  effective  sand  size,  bulk  density, 
filter  drag  coefficient,  single  grain  collection  efficiency,  and 
solidarity  factor  are  determined  for  each  bed  depth  increment  as  shown 
in  Figure  1.  This  is  done  to  account  for  the  gradation  occurring  when 
a sand  filter  is  backwashed.  The  result  is  effectively  constructing  a 
ten-layer  sand  filter  of  which  each  layer  is  ten  percent  of  the  over- 
all sand  bed.  The  effective  sand  size  for  each  bed  increment  is 
determined: 


d (I)  = exp[ln  d + f(I)  In  °q] 

where  d^s(I)  is  the  effective  sand  grain  size  for  bed  increment  I 
(cm),  dgs  is  the  average  effective  sand  grain  size  for  the  filter  (cm), 

°g  is  the  log-normal  standard  deviation  for  the  sand  grain  size  distri- 
bution (cm),  and  f(I)  is  the  log-normal  statistical  coefficient.  This 
coefficient  has  values  corresponding  to  the  five  through  ninety-five 
percent  size,  in  a log-normal  distribution,  that  permit  calculation  of 
the  sand  grain  size  that  corresponds  to  the  respective  percentage 
using  the  general  form  of  equation  7. 

Because  sand  that  is  distributed  log-normally  in  size  will 
settle  according  to  Stokes'  law,  sand  filters  exhibit  a size  gradient 
from  top  to  bottom.  This  gradient  is  created  when  the  sand  filter  is 
backwashed  and  allowed  to  settle.  The  large  sand  grains  have  a 
greater  terminal  velocity  than  do  the  smaller  grains.  For  this  reason 
sand  grains  in  the  filter  bed  increase  in  size  from  top  to  bottom. 

Since  particle  bed  density  increases  as  particle  size  decreases,  the 
bulk  density  of  the  bed  decreases  from  top  to  bottom.  To  account  for 
this  increase  in  bulk  density,  a linear  approximation  was  used  in  the 
model  between  the  bulk  density  of  the  ten  percent  size  and  the 
average  bed  bulk  density.  These  two  parameters  were  inputs  to  the 
model;  the  linear  approximation  calculates  the  gradient  between  sand 
layers  and  assigns  the  appropriate  bulk  density  to  each  layer  of  the  bed. 


mm  ■ ’•••  - •• 


The  filter  dray  term  for  each  sand  layer  is  calculated: 
1.84  x IQ'3  P(I) 


K2s(I)  = 


dp<;(I)2  P 2 [1  - (pRc(I)/Pc)]: 


where  K?S(I)  is  the  drag  term  for  sand  layer  I (sec/cm),  PBS  (I) 
is  the  Bulk  density  of  sand  layer  I calculated  from  the  above 
described  linear  approximation  (g/cm3),  des(I)  is  the  effective 
sand  size  in  sand  layer  I calculated  by  equation  7 (cm),  and  ps 
is  the  density  of  a discrete  sand  grain  (g/cm3).  The  density 
of  the  discrete  sand  grains  is  constant  for  single  media  filters. 

The  single  grain  collection  efficiency  is  determined  from  the 
modified  Friedlander  (References  12,  13)  equation;  however,  for 
the  purposes  of  this  study,  the  rate  of  capture  was  assumed  to  be 
unity.  The  single  grain  collection  efficiency  expression  may  be 
written: 


n (I) 
ICD<|>  ' 
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where  njflud)  is  the  modified  Friedlander  (References  12,  13) 
collection^efficiency  for  a single  sand  grain  in  bed  increment  I 
and<)>is  rate  of  capture  for  particles  colliding  with  a sand  grain, 
which  is  assumed  to  equal  unity. 

The  last  calculation  made  for  the  layered  filter  is  the  soli- 
darity factor.  This  is  perhaps  the  most  influential  parameter  in 
the  graded  media  filter  and  has  been  the  basis  for  the  develop- 
ment of  reverse  graded  filters  in  recent  years  (References  6,  14- 
18,  20  - 24).  As  discussed  in  reference  to  bulk  density,  a sand 
filter  media  will  settle  after  fluidization  with  a sand  grain  size 
gradient  from  the  top  to  the  bottom  of  the  bed.  The  grain  size 
increases  along  this  gradient.  This  is  unfortunate  because  the 
collection  efficiency  of  a single  sand  grain  is  inversely  propor- 
tional to  its  size;  hence  the  top  layer  of  sand  collects  suspended 
particles  most  efficiently. 


There  are  two  ramifications  from  this  phenomena:  (1)  because  the 
high  efficiency  grains  make  up  the  upper  layers  of  the  media, 
particles  that  penetrate  the  top  layer  are  less  likely  to  be  col- 
lected in  each  successive  layer  and  (2)  the  rapid  accumulation  of 
the  greatest  portion  of  the  suspended  particles  on  the  upper  sand 
layer  increases  the  rate  of  suspended  particle  formation  and  hence 
the  rate  of  increase  in  pressure  drop  with  time.  Under  high  head- 
loss  situations,  the  surface  of  the  cake  may  collapse  or  the  par- 
ticles, if  flocculant  in  nature,  may  be  crushed.  If  this  occurs, 
floe  break-through  may  result  since  the  low  efficiency  lower  sand 
layers  may  be  unable  to  collect  the  small  particles  that  have 
escaped  from  the  upper  sand  layers  and/or  the  suspended  particle 
cake. 

The  concept  of  the  reverse  graded  filter  allows  the  particle 
size  of  the  graded  media  to  decrease  from  the  top  of  the  filter 
bed.  In  this  instance  greater  media  penetration  can  be  obtained. 
Thus  the  rate  of  pressure  buildup  with  respect  to  time,  for  a 
given  system,  should  be  reduced. 

The  above  discussion  emphasizes  that  the  log-normal  standard 
deviation  of  sand  size  in  a sand  filter  is  a very  influential 
parameter.  The  uniformity  coefficient  has  long  been  used  in 
description  of  filter  sand;  this  term  is  a rough  approximation  of 
the  log-normal  standard  deviation.  The  solidarity  factor  for  each 
sand  layer  is  determined  from  the  characteristics  of  that  layer 
according  to  the  solidarity  equation  presented  previously.  This 
expression  may  be  written: 

Sf(I)  . °-'90"  Ls  PBS(1> 

ps  ^es  (I)  (10) 

where  S^(I)  is  the  solidarity  factor  for  sand  layer  I (unitless) 
and  is  the  overall  sand  bed  depth  (cm).  As  shown  in  Figure  1, 
equations  7 through  10  are  repeated  for  each  sand  layer  from  the 
top  to  the  bottom  of  the  sand  filter. 

After  the  characteristics  of  the  clean  sand  filter  are  estab- 
lished, the  characteristics  of  the  suspended  particles  must  be 
determined.  As  shown  in  Figure  1 the  first  of  these  characteristics 
to  be  determined  is  the  Friedlander  equation: 


nICD*  a 6 [ _£ 


u depC9 


u d p 
ep  g 


where  icD<t>  is  the  modified  Friedlander  (References  12,  13)  collection 
efficiency,  de  is  Chen's  (Reference  10)  effective  size  as  a collector 
(this  is  calculated  using  equation  4,  using  the  surface  mean  suspended 
particle  size  and  the  arithmetic  mean  suspended  particle  size;  the 
dimension  is  cm),  dnp  is  the  arithmetic  suspended  particle  size  (cm), 
and  <t>'  is  the  probability  of  a successful  attachment  upon  collision 
which,  for  the  purposes  of  this  study,  is  assumed  to  be  unity.  At 
this  point  in  the  sequence,  it  is  necessary  to  make  an  assumption 
with  respect  to  the  sand/suspended  particle  interface  at  the  top  of 
the  filter.  Since  the  top  grains  of  sand  are  the  smallest  in  the 
size  distribution,  they  are  assumed  to  supply  the  nucleus  for  sus- 
pended particle  cake  formation.  Thus,  it  is  assumed  that  the  surface 
of  the  filter  is  composed  of  a one-particle-deep  layer  of  suspended 
solids;  the  diameter  of  the  particles  is  equal  to  the  calculated 
equivalent  suspended  particle  diameter  when  acting  as  a collector. 
Thus,  the  initial  cake  solidarity  factor  is  calculated  using  the 
effective  suspended  particle  diameter: 

, 1.9099  Lr  PRp 

Sf  ■ j— C-  BP  (12) 

Pp  dep 

• 

where  Sf  is  the  initial  solidarity  factor  of  the  suspended  particle 
cake  (unitless),  Lc  is  the  initial  cake  thickness  equal  to  d«n  (cm), 
the  effective  particle  diameter  as  a collector  according  to  {Reference 
10)  equation  (cm),  p^  is  the  bulk  density  of  the  suspended  particle 
cake  (g/cm3),  and  p_  is  the  density  of  a discrete  suspended  particle 
(g/cm3).  p 

Once  the  calculations  in  equations  1 through  12  are  accomplished, 
numerical  integration  of  pressure,  flow,  and  efficiency  as  a function 
of  time  can  be  initiated.  This  is  done  by  incrementing  time  so  that 
each  increment  is  insignificant  with  respect  to  the  final  accuracy 
required  as  shown  in  Figure  1. 
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Variables  that  are  affected  by  the  changing  characteristics  of 
the  filter  media,  as  it  accumulates  suspended  particles  on  its 
surface  and  in  the  sand  itself,  are  adjusted  to  reflect  this 
change  upon  flow,  time,  pressure,  and  efficiency  during  each 
iteration.  These  modifications  are  accomplished  sequentially  as 
the  graphical  integration  procedure  increments  time. 

In  order  to  facilitate  understanding  of  the  integration 
incremental  sequence,  the  reader  should  refer  to  Figure  1 as 
its  mathematical  sequence  is  described.  Suspended  particle 
cake  thickness  is  incremented: 


CTutTnr(t-tT) 

-c(t)  = Lc  ( t-tj ) + J-.1  C L 


(13) 


BP 


where  L~(t)  is  suspended  particle  cake  thickness  (cm)  at  time  t 
(sec),  Lc(t-tj)  is  the  suspended  particle  cake  thickness  at  the 
previous  time  increment  ( t j is  the  incrementing  unit),  nc(t-tr) 
is  the  collection  efficiency  for  the  suspended  particle  cake  during 
the  previous  time  increment,  and  u is  fluid  velocity  (cm/sec). 

Once  the  new  cake  thickness  is  determined,  aP  (t)  pressure  drop 
caused  by  the  cake  (formed  from  the  suspendedparticles)  can  be 
determined: 


3 


\ 


I 


| 


:: 


APc(t)  = K2p  U p I (t) 

Bl  C (14) 

where  APc(t)  is  the  pressure  drop  caused  by  the  suspended  particle 
cake  (g/cm2)  at  time  t,  K2p  is  the  suspended  particle  unit  depth 
drag  coefficient  (sec/cm) ,Hand  Lc(t)  is  the  suspended  particle  cake 
thickness  at  time  t. 

In  addition  to  the  pressure  drop,  the  new  solidarity  factor  for 
the  suspended  particle  cake  can  be  determined  at  this  point  as 
follows: 


$.  = '-9099  pBCLc<*> 

D (T 
Pp  °ep 


(15) 


I 

where  Sf(t)  is  the  new  suspended  particle  cake  solidarity  factor  for 
the  increased  cake  thickness  (dimensionless)  Lc(t ) at  time  t.  At 
this  point  the  collection  efficiency  for  the  suspended  particle  cake 
can  be  calculated: 
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-sf  (t)n 


(16) 


n (t ) = 1 - exp 
%0 


ICO<t> 


where  \(t)  Is  the  suspended  particle  cake  collection  efficiency  at 
time  t and  nl.j.  is  the  modified  Friedlander  (Reference  12,  13) 
suspended  particle  upon  suspended  particle  collection  efficiency. 

Once  the  f'luid/particle  suspension  has  passed  through  the  particle 
cake  on  the  surface  of  the  sand  media,  the  concentration  of  suspended 
particles  in  the  fluid  is  reduced  by  the  percent  collected.  For  this 
reason  the  suspended  particle  concentration  reaching  the  top  of  the 
sand  bed  must  be  calculated: 

CI2  = Cj[l  - nc(t)]  (17) 

where  CI2  is  the  suspended  solids  concentration  reaching  the  sand 
surface  and  Cj  is  the  initial  suspended  solids  concentration. 

At  this  point  another  sequence  is  necessary  to  calculate  the  con- 
centration gradient  of  the  suspended  solids  in  the  fluid,  as  the  fluid 
passes  through  the  sand  bed.  This  is  accomplished  by  determining  the 
mass  of  solids  retained  in  each  sand  layer  as  a function  of  time. 

Flow,  time,  pressure,  and  efficiency  inter-relationships  are  reeval- 
uated for  each  successive  time  increment  for  all  sand  layers.  The 
first  term  calculated  in  this  sequence  is  the  sand  layer  efficiency: 

VD  - 1 - exp  [sf(»nico/lT|  (18) 

where  ns(I)  is  the  efficiency,  Sf(I)  is  the  solidarity  factor,  and 
nsU)  is  the  efficiency,  Sf(I)  i$  the  modified  Friedlander  (11,  12) 
single  grain  collection  efficiency;  all  of  these  are  for  the  I™ 
sand  level . 


The  next  calculation  predicts  the  pressure  drop  across  the  Ith 
sand  level: 


APS(I)  = 


K2S^)  * Ls  W1* 

10. 


(19) 


where  aP  (I)  is  the  pressure  drop,  L is  the  total  sand  bed  depth, 
pdc(I)  is  the  sand  bulk  density  for  the  ith  sand  level,  and  u is 
fTTtration  velocity. 
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The  next  calculation  determines  the  new  bulk  density  of  the  sand 
level  after  suspended  particles  have  been  collected  during  the 
integration  increment.  This  calculation  may  be  expressed: 


PBsUJti  = PBS(I)  + CI2(l)u  tt  nICDd)(I) 
(Ls/10. ) 


(20) 


where  PRs(*)ti  is  the  new  bulk  density  for  the  next  time  increment, 
tj  is  tne  time  increment,  and  nicpd>U)  is  the  efficiency  of  the  sand 
layer  during  the  current  increment  as  calculated  by  equation  18. 


Once  the  new  sand  bed  characteristics  of  the  layer  are  determined, 
the  filter  drag  and  solidarity  factor  are  modified  to  consider  the 
accumulated  suspended  particles.  These  modifications  can  be  accom- 
plished once  the  value  of  the  new  bulk  density,  determined  by 
equation  20,  is  calculated.  The  other  variables  affecting  the  drag 
coefficient  and  the  solidarity  are  constant  within  the  integration 
operation.  Filter  drag  is  determined  by  taking  the  previous  inte- 
gration increment  drag  term  and  adding  the  drag  increase  contributed 
by  the  newly  entrapped  suspended  particles.  This  is  calculated: 

«*,(!)«  - K2s(.)»BS(»  * (n) 

pBS(I)ti 


where  K2$(I)ti  is  the  new  filter  drag  term  for  sand  layer  I. 

The  new  solidarity  factor  for  the  sand  layer  is  determined 
directly: 

= 1-9099Ls  PBS(I)ti  (J 

10.  ps  des ( I ) 

where  Sf ( I )ti  is  the  new  solidarity  factor  for  sand  layer  I. 

Once  these  modifications  are  made  to  the  1^  layer  sand  bed 
characteristics,  the  values  of  interest  are  accumulated.  These 
values  include  Ith  level  efficiency,  pressure  drop,  and  efficiency 
from  the  top  of  the  sand  bed  through  the  level.  The  efficiency 
parameter  is  determined  by  the  accumulation  of  efficiency  as  the 
fluid  flows  through  successive  layers  of  sand.  It  should  be 
recognized  that  the  accumulation  of  efficiency  terms  must  express 
the  successive  increase  as  a function  of  the  initial  water  quality, 
not  as  a function  of  the  water  quality  entering  the  layer  under 
consideration.  For  this  reason  the  general  addition  efficiency 
expression  is  written: 
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ns  = 1 - 0.  - nSA)  (1.  - n j ) (23) 

where  n$  is  the_total  efficiency  from  the  top  of  the  filter  through 
the  Ith  level,  n™  is  the  total  efficiency  through  the  level 
immediately  above  the  ith  level,  and  m is  the  efficiency  in  the  Ith 
level.  1 

Unlike  efficiency,  the  total  pressure  drop  at  any  point  in  the 
filter  is  additive.  Thus,  the  total  pressure  loss  can  be  expressed: 

aPs  = aPsa  + aPi  (24) 

where  aPs  is  the  total  pressure  at  the  bottom  of  the  Ith  level, 
aPSA  the  total  pressure  drop  through  the  level  immediately  above 
the  Ith  level,  and  aPj  is  the  pressure  drop  in  the  ith  level.  The 
last  parameter  that  requires  determination  in  this  incremental 
evaluation  of  the  sand  bed  is  the  reduction  of  the  suspended  particle 
concentration  through  the  sand  layers  above  and  including  the  Ith 
sand  layer.  This  is  calculated: 

°I2(I)  ' CI2  ' "s’  (25) 

where  Cjod)  is  the  new  suspended  solids  loading  in  the  fluid  after 
the  fluid  has  passed  through  layer  I and  C j2  is  the  suspended  particle 
loading  at  the  top  of  the  sand  bed  after  it  has  passed  through  the 
suspended  particle  filter  cake. 

Equations  18  through  25  are  repeated  through  each  of  the  ten  sand 
layers  and  the  values  of  all  parameters  are  accumulated  as  prescribed 
by  their  respective  equation.  Once  this  is  accomplished,  as  indicated 
by  Figure  1,  the  overall  pressure  drop  and  efficiency  for  the  time 
increment  is  determined,  and  the  average  filtration  cycle  efficiency 
is  determined.  The  overall  pressure  drop  is  calculated: 

aP^  = aPc  + aPs  (26) 

where  APt  is  the  total  pressure  drop  across  the  filter,  aPc  is  the 
pressure  drop  across  the  suspended  particle  cake  as  calculated  by 
equation  14,  and  Ps  is  the  overall  pressure  drop  across  the  sand  bed 
as  determined  by  equation  24.  The  incremental  efficiency  at  any 
time  t is  determined: 

njt)  = 1.  - [1.  - nc ( t ) ] [1.  - ns(t)]  (27) 

where  n*(t)  is  the  total  filter  efficiency  at  time  t into  the  cycle, 
nc(t)  and  ns(t)  are  the  efficiencies  of  the  suspended  particle  cake 
and  the  sand  bed  during  the  last  time  increment,  respectively. 
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Both  of  these  terms,  therefore,  represent  their  respective  efficien- 
cies at  time  t into  the  filtration  cycle. 


In  addition  to  the  level  of  filter  efficiency  at  time  t,  the 
average  efficiency  from  the  start  of  the  filtration  cycle  to  time  t 
is  of  interest.  This  variable  is  calculated  by  accumulating  the 
incremental  efficiencies  and  then  dividing  by  the  number  of  increments. 
It  is  apparent  that,  unless  floe  breakthrough  occurs,  the  filter 
efficiency  will  increase  as  a function  of  time.  For  this  reason 
the  average  cycle  efficiency  will  always  be  less  than  the  current 
efficiency  at  time  t. 


Equation  sequence  1 through  27  is  repeated  for  the  time  increments 
until  the  desired  time  (or  pressure  drop)  for  filter  backwash  is 
reached.  This  time  may  well  be  determined  by  a trial -and-error  method 
to  facilitate  filter  operation  between  a residual  pressure  drop  and 
some  maximum  desired  pressure  drop  that  allows  the  most  economical 
operation  of  the  filtration  system. 

It  would  be  extremely  difficult  to  perform  these  manipulations 
manually.  This  sequence  is  complex  in  that  it  considers  many  filter 
design  parameters  that  are  not  normally  taken  into  account  in  more 
classical  design  sequences.  Because  use  of  the  digital  computer  has 
become  widespread,  application  of  complex  concepts  and  lengthy  design 
sequences,  previously  unused  because  of  their  arduous  nature,  can  be 
easily  accomplished;  such  is  the  case  with  this  design  technique  for 
sand  filtration  systems.  Since  the  design  engineer  can  interact  with 
the  computer  in  the  time-sharing  mode  he  can  consider  a multitude  of 
different  concepts  that  were  previously  too  complex  or  too  time 
consuming  for  evaluation. 
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